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Abstract. The two body problem in a scalar theory of gravity is investigated. We focus on the closest theory to General Rel- 
ativity (GR), namely Nordstrom's theory of gravity (1913). The gravitational field can be exactly solved for any configuration 
of point-particles. We then derive the exact equations of motion of two inspiraling bodies including the exact self-forces terms. 
We prove that there is no innermost circular orbit (ICO) in the exact theory whereas we find (order-dependent) ICOs if post- 
Newtonian (PN) truncations are used. We construct a solution of the two body problem in an iterative (non-PN) way, which 
can be viewed as a series in powers of (v/c)^. Besides this rapid convergence, each order also provides non-perturbative in- 
formation. Starting from a circular Newtonian-like orbit, the first iteration already yields the 4.5 PN radiation reaction. These 
results not only shed light on some non-perturbative effects of relativistic gravity, but may also be useful to test numerical 
codes. 
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INTRODUCTION 

The two-body problem in GR is both an important and difficult issue. Various approach have been investigated. One 
generically needs some assumptions and simulations (numerical relativity) or perturbatives techniques (PN expansion, 
e-g El)- In both cases it remains difficult to quantify the errors one makes. This point is however crucial in order to 
interpret the coming-soon data. 

It should thus be very interesting to test the numerical codes and the assumptions usually made in the two-body 
problem within a toy model of GR. One may also hope that some generic non-perturbative effects may be found, that 
could also occur in GR. This toy model should be simple enough to treat analytically (at least partially) the two body 
problem, but should also be as close as possible to GR. Metric scalar theories of gravity are perfect candidates of such 
toy models. 

This idea is not a new one. It has already been used in order to test numerical codes, see e.g or, more recently, 
the validity of the Quasi Equilibrium (QE) scheme 101 ■ However we wish here to focus on a particular theory of 
this class (namely Nordstrom's one), which is, to our mind, the best motivated toy-model of GR, and which has not 
been used in the literature. An important point is that Nordstrom's theory respects the Strong Equivalence Principle 
(SEP) and thus enables us to consider constant-mass points-particles, whereas this cannot be justified in others models 
considered in the literature. We also want to go further in the analytical resolution of the two-body problem, by using 
one of the great advantages of Nordstrom's theory: the linearity of the field equation. 

In Section II we recall the basic features of such metric scalar theories and notably the particular case of Nordstrom's 
theory. In Section III, we solve exactly the metric in terms of arbitrary massive point particles. We then derive the 
exact equations of motion of such bodies which notably involve finite self-force terms. 

In Section IV, we restrict to the equal-mass two-body case and investigate the circular motion. We derive the 
analytical relation between the radius and the orbital velocity. We compute the exact energy of the binary in such 
a circular configuration and show that the theory does not exhibit an innermost circular orbit (ICO), whereas such 
ICOs are found if one works with PN truncations of the exact energy. 

In Section V, starting from this circular solution, we show that the inspiral motion can be found order by order 
in power of (v/c)^, where v is the velocity of the bodies. This behavior is directly related to the leading multipolar 
emission of the system, namely the quadrupolar one. Furthermore, this expansion is not a PN one since each orders 
yield non-perturbative information on the motion. 

In Section VI, we show some analytical and numerical results. We give the exact rate of change of the velocity 
up to 4.5 PN and the radius up to VPN. We also plot these two quantities and by numerical integration, we plot the 



corresponding inspiral motion. This paper is a short version of the self-contained article Jit] in which more technical 
details can be found. 



A QUICK LOOK AT NORDSTROM'S THEORY OF GRAVITY 

We first introduce briefly metric scalar theories of gravity. The action reads 

=^ = -^-7; / \/^ T7''^5^^av<P + J5^matter[matter;|^v=A2((p)T7^v], (1) 
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where A{(p) is a given function of the scalar field, characterizing the theory. The action of matter, ^matter' ^ 
functional of all matter fields, assumed to be minimally coupled to the "physical" metric = ^^(9')^7/iv- These 
theories thus respect the weak equivalence principle. Throughout this paper, we will use the sign conventions of |6], 
and in particular the mostly-plus signature — h ++. Throughout the paper we will use coordinates for which the flat 
metric T7|jv takes its fundamental form. We use bold-faced symbols to denote the three vectors of the Minkowskian 
geometry. 

Since the physical metric is conformally related to the flat one, there is no coupling of the photon to the scalar field. 
In any theory of the type there is thus strictly no light deflection, and all of them are ruled out by experiment. They 
anyway share many feature with general scalar-tensor theories J?!] and one of them even satisfies the SEP. 

The field equations deriving from action ^ read 

AnG 1 

^^^t<P = ~—AH<p)A'{<p)T (2) 

V^,T^^=Q, (3) 

where t^^ is the physical stress-energy tensor, f = gfivf^^ its trace, denotes the covariant derivative with respect 
to (w.r.t) the physical metric g^y, and Dflat d' Alembertian operator. The scalar curvature may be written as 

R^'-^A'H,)T-6^^n,d.,. (4) 

Therefore, if and only if A((p) = (p, the theory admits a purely geometrical and generally covariant formulation 

R^^^f,cf,yp=0,V^T^''^0 (5) 

where C denotes the Weyl tensor. Using the same reasoning as in J^, this geometrical formulation suffices to prove 
that the SEP is satisfied. Among all others theories of type OJ, this particular one characterized by A{(p) = (p is thus 
the closest one to GR. In the following, it will be referred to as Nordstrom's theory of gravity ifToll (see also ifTlll and 
the review 1 12]). 

To our knowledge, there exist only two gravity theories satisfying the SEP 10]: GR and Nordstrom's theory of 
gravity. The SEP notably means that the gravitational binding energy of a body contributes the same to its inertial 
and gravitational masses, so that strongly self-gravitating bodies fall in the same way as test masses in an external 
gravitational field (up to self-force effects which can play a major role in the dynamics). This fact therefore allows us 
to describe massive bodies as constant-mass point particles without worrying about their internal structure. In all other 
theories of type ([T]i, with a non-linear matter-scalar coupling function A (9), violations of the universality of free-fall 
of self gravitating bodies would appear already at the first PN level. 

In the following, we will be interested in matter consisting of massive point particles of coordinates z^, whose 
Lagrangian read L — — Y.A'^AcA{(p)ds where ds^ = rj^ydx^dx^ is the line element of the flat metric. Then the field 
equation reads : 

Oflat'PW = ^A'((p)£mA / (x^ -z^ (t^)) ds, (6) 

C A •' 

where T is the proper time along worldlines. Here we notice that Nordstrom's theory, corresponding to A(^) = ^, gives 
a linear equation, so that (p is simply the sum of a constant and of the separate contributions of each point particle. 
This a great theoretical advantage over GR. The linearity of the field equation together with the fact that the SEP holds 
in Nordstrom's theory motivates our study of this particular theory. 



DYNAMICS OF N POINT-PARTICLES; THE TWO BODY PROBLEM 



Let us consider point particles whose motion is assumed to be known. We parameterize their positions z^(t^) by 
their proper time Ta- Then the linear field equation Eq. (|6]l is solved with 
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where 

pf\x) = ^4[x'^){x'-zl{T'A)), (8) 

is a scalar distance between x and za('!a)' where u is the four-velocity and is the retarded proper time given by the 
intersection of the A-th particle world line and the past light cone of x. Following I I3i1 . we refer to it as the retarded 
distance. Let us stress that, contrary to GR or to any theory of type ([TJ with A ^ Id, the gravitational field is thus 
exactly known. 

On the other hand, we now look for the equation of motion of point particles in a given external field. Equation (O 
shows that test masses are following geodesies of the physical metric. However, in the two comparable mass body 
problem, one cannot neglect the effect of the proper field on the motion (the so-called self-force). Actually, Eq. (|7]l 
shows that the self-field is singular on bodies world lines, so that their exact equation of motion (including the self- 
force) is not yet defined. 

The same problem occurs in classical electromagnetism. Dirac addressed this issue in EM by using the local 
conservation of energy-momentum in a small 3-tube surrounding the world lines, and derived the correct (known as 
the Lorentz-Dirac) equation of motion. In |5] we adapt this procedure to our specific case using Eq. Q, following lITsil 
and ifisll . The equation of motion of body A in a smooth external field (p^^^ finally reads 

^ ((l + (p^^^) Mv) + ^ (li^Mv - iiv) + <3v(p^''^ - 0, (9) 

where is the mass of body A, u denotes its four- velocity and the dot means the derivative w.r.t the proper time. As 
far as the A^-body problem is concerned, one should reintroduce the indices and write the external field with the help of 
Eq. (|7]i. The self-forces term in Eq. (|9|l are those which are proportional to niA (because neglecting the self-force means 
neglecting the mass of the particle), and we see that they are third derivative of the position, as in the Lorentz-Dirac 
equation. 

Let us specialize these results to the two body problem. The above manifestly covariant form of the equation 
of motion (EOM) can be written in a more useful way in terms of cartesian coordinates {ct,x). We write z^(f) = 
{ct,ZA{t)), and similarly for B. Let VA{t) — dzA/dt be the three-velocity. We define its velocity as J3^(f) = VA{t)/c, 
and the Lorentz factor by (1 — /S^)^'/^. For any position of body A at time f, there is a unique retarded position 
of body B given by the intersection of the past light cone of z^(f) and the world line of B. We shall write it as 
Zg(fi-et) = (cfretiZB(fret))' where fj-gt is the retarded time. The retarded distance is then 

pf\zA{t)] - -«^(fret) {zl{t) -z^(fret)) , (10) 

and is simply denoted p^^- With the above notations, it is straightforward to show that the dynamics of body A at time 
t is given by : 



Yl{\-GmB/pf'c^) 
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(11) 



where all quantities have to be taken at time t. Now the dot means the derivative w.rt f/me, d/dt. The gradient of p^^^ 
is the one w.r.t the position of body A (see ^ for useful formulaes). The left hand side is just the acceleration of body 
A, and we recognize in the first term of the right hand side a Newtonian-like term, responsible of a force proportional 
to mAMB, whereas the second one is the self-force term, proportional to m^. The exact equation of motion for the two 
body problem is thus known. 

In the following, we shall restrict to the case of equal-mass bodies. The EOM then trivially comes from Eq. ( fTTI ) 
with niA = ms = ni. In that case we have a symmetry that ensures the existence of a Lorentz frame in which the center 
of mass is always at rest. In this frame, the three-position vectors w.rt the center of mass obey the law 7^(0 = ^ 



for all time f. Furthermore the motion can be shown to be plane (see [5]). As a consequence, the motion can be 
characterized by only two variables whose choice is free. In the following we choose r[j3 (f)] and j3 (f), where j3 (f) is 
the norm of the velocity of A (or B) and r is the norm of za (or B). The reason of this choice is that the gravitational 
radiation (and therefore the non-trivial temporal evolution of the binary system) is encoded in the variations of j3(r). 
The angular velocity w is easily expressed in terms of these variables. 

It must be stressed that our dynamical equation is not a complete set of equations. Indeed, we have to write the 
kinematical equations that determine the retarded position of B. If we define dt^^^ = t — fj-gt W^^^ the oriented, 
positive, retarded angle between Zb (fret) ™d — z^(f ), where the positive orientation is chosen to be the one defined by 
the angular velocity vector, these (implicit) equations are 

c8f^' = ^ P- [j3 (f )] + 2r[^ (f)]r[j3 (fret)] cos w'^' + [H (fret)] , w'^' = J' (o{t')dt' . (12) 

This set of four equations Eq. (fTTl l (with niA — ins ~ m) and Eq. (fT2l l. which we will refer to as the equations of motion 
(EOM), are now sufficient to determine the entire motion. In the next section, we look for a circular motion. This 
unphysical motion will be used as an initialization of an iterative method that will construct order by order the actual 
inspiral motion, see Section V. 



THE CIRCULAR CONFIGURATION, ITS ENERGY AND THE ICO 

We now restrict ourselves to equal-mass binaries. In order to obtain a stationary, circular solution, we have to use the 
time symmetric Green function when solving the field equation, i.e., we have to consider as much outgoing waves as 
ingoing ones. We show in jsll that the circular solution then reads 



Gm 1+213^- jg^cos y/Q - 2 (j3 + jg^) sin ' ^ 
^ 47j32cos(f ) (l+)3sin(f )) 



where the radius ro denotes the distance of one body to the center of mass, and j3 —v/cis the constant orbital velocity 
of the bodies. In the literature the notion of separation is more often used and is given here by 2rQ. Here m denotes the 
mass of each body, and the retarded angle \j/o is given by the implicit equation 



f^^^cos^^J. (14) 

The index "0" simply refers to the circular case. Note that to the leading order, the radius reads ~ rNewton = 
Gm/A^^c^ just as in Newton's theory. In order not to show this leading but trivial part of the function tq, we plot 
in the right panel of Fig. ^ the ratio ro lri^i,^^^on as a function of j3 . We see that it goes to zero if the velocity goes to 1 . 

Since Nordstrom's theory shares many features of the general relativity, it is interesting to look for the exact 
expression of the energy of the binary being in a ckcular motion, in order to investigate in a non-perturbative way 
the existence or not of the Innermost Circular Orbit (ICO), which is an important notion in numerical relativity even if 
it does not correspond to any physical observable. 

In Jst] we compute the energy of the binary using the Fokker action associated to the initial action of the theory. The 
result is surprisingly simple since it reads : 

^(^) = ^, (15) 

which have also been found independently by Friedman and Uryw ifTsIl in a related context. The energy of the binary 
thus goes to zero in the ultra relativistic limit, so meaning that all the initial energy of the binary has been carried away 
by radiation. Since the limit v = c is also the limit ro = 0, we shall consider that the two bodies have melt each other 
This remaining point particle has a vanishing energy, and thus does not actually exist. The two initial bodies have 
thus been entirely evaporated into gravitational radiation. It must however be stressed that this conclusion would hold 
only if the physical motion were a succession of circular orbits of decreasing radius, which is obviously not the case. 
Actually the inspiral motion can be seen as such a succession only if the damping timescale is much greater than the 
orbital one, an assumption that breaks down in the relativistic regime, as we will show in Section VI. 

The ICO is defined as the separation of the companions for which the energy is minimum, if such a point exists 
IitIi . Now, if the separation decreases and becomes smaller than the ICO, the motion cannot remain (quasi) circular 
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FIGURE 1. The PN ICO's up to 38 PN. ICO's do not necessarily exist at each PN order. The series seem to converge whereas 
this can not be the case. ICO's found in GR are also shown (extracted from I1I6I ). 



unless energy of binary grows, which is impossible because of gravitational radiation. Passing through the ICO may 
therefore be the definition of the beginning of the plunge phase. Since 
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does not vanish (unless j3 = 0, that is r = 00), there is no ICO in Nordstrom's theory. However, if one now works in a 
PN truncation, one gets e.g. at the second order 
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This 2PN energy has a turning point for a separation being roughly l.OSGm/c^. It thus shows that order-dependent 
ICOs can be found although the non-perturbative ICO does not exist. We find such PN ICOs at 2, 3, 5, 6, 8, 9,. . . PN 
orders. In Fig. ([T]i we have plotted up to very high PN orders (38 PN) the position of the ICO and placed the points 
obtained in general relativity (extracted from 01611 ). In this figure, we plot the angular frequency of the binary at the 
ICO as a function of the energy of the binary at the ICO, normalized to the total energy : [E^qq — 2mc^) / [Imc^) (see 
S). 

The ICOs we find here are quite faraway from those found in GR (by numerical study or PN expansions). The most 
striking feature of this plot is that the PN ICOs seem to converge, whereas there is no ICOs in the exact theory. 



CONSTRUCTION AND CONVERGENCE OF A PERTURBATIVE SOLUTION 

Beyond the circular motion, which is not physical unless advanced and retarded waves are considered, we shall now 
construct an inspiral motion, solution of the only retarded EOM. When deriving a perturbative solution to the EOM, 
one find convenient to look for the orbital acceleration /3 as a function of j3 

m=fm% (18) 

which, together with the radius as a function of )3, completely characterizes the motion. The perturbative expansion of 
the EOM goes as follow. We look for function / and r in power of a parameter e, as / = /o + e/i + . . . + e"/„ + . . . 
and similarly for r, where the ft and r,- are unknown function of /3, and e is an arbitrary parameter. We choose /o = 
in order to recover an expansion whose 0-th order corresponds to the circular motion. It is important to note that this 
choice is analytically unmotivated, but rather corresponds to an initial condition assumption. 



We have proven in jst] that this expansion converges in powers of j3^, which means that the property fi ^ and 
r,- ^ is satisfied when e = j3^. Moreover, this amplitude is directly related to the first non-vanishing multipolar 
radiation of the binary, namely the quadrupolar one : E/E ^ j3^, where E is the energy of the binary. Note also that 
this expansion is not a PN one since each functions /, and r, are found to possess a complicated series expansion, and 
are not just monomials. For instance the 0-th order provides an analytical formula linking the radius of the circular 
motion to its orbital velocity, see Eq. ( fT3] l. This method thus gives to each order many non-perturbative information. 



ANALYTICAL AND NUMERICAL RESULTS 

We have run the previous algorithm up to the second order It means that we have derived the exact, analytical 
expression of /i, /2, ri and r2. We do not write explicitly these functions here since their expression are quite heavy, 
but a simple code enables to find them. These functions / and r fully characterize the motion, so that we will refer to 
these first and second order as the first and second order motion. 

We can write the explicit and exact PN expansion of these two solutions. We show in ^ that the first order of 
perturbation give both the acceleration and the radius exactly up to 4.5 PN, whereas the second order provides a 
correct radius up to 7 PN. These results are : 
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Let us now focus on the entire non perturbative functions. In the left panel of Fig. we plot the orbital acceleration 
as a function of j3 . We show the behavior of this acceleration for the first and second order motion, and we also plot 
the 4.5 PN expansion of it given in Eq. ( fT9] l. We also plot in the right panel of Fig. (|2]l the radius as a function of /3, 
normalized to the Newtonian radius ^^^ewton = Gm/Afi^c^. We show the circular radius, the first and second order 
one, and the 7PN radius given by Eq. ( |20l ). 
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FIGURE 2. (left) The orbital acceleration as a function of /3: the first order and second order solutions, compared to the 4.5 PN 
one. (right) The radius of the orbit as a function of j8. We plot the circular solution, and the first and second order corrections to it. 
The radius at 7 PN is also given. 

In Jsl], we also investigate the so-called Quasi-Equilibrium (QE) scheme which has been extensively used in 
GR (see, e.g, lfT9ll . ll20ll ) and discuss its validity in Nordstrom's theory. This approximation holds if the dynamical 



(damping) timescale is much greater than the orbital period, so that the gravitational radiation can be neglected at first 
approximation, and the motion is quasi-circular. The left panel of Fig. ^ represents the ratio (at first order) of the 
orbital period over the dynamical damping timescale. We thus see that the QE approximation is an excellent one in 
the non-relativistic and in the ultra-relativistic regime. However, in the intermediate regime ()3 ^ 0.6) the dynamical 
timescale becomes comparable and even shorter than the orbital period, so that the orbit is highly non-circular and the 
approximation underlying the QE-scheme is strongly broken. This is the plunge phase of the binary. 

A numerical integration of Eq. ( fTSl l yields the behavior of the norm of the velocity j3 as a function of time, plotted 
in right panel of Fig. (O. The three curves corresponds to the 4.5 PN result, the first order and second order result, 
starting at )3 =0.1. The three results are very similar until /3 ~ 0.4. A major difference, however, between the PN 
solution and the non-perturbative solutions is that it cannot take into account the pole of special relativity, so that 
the PN solution diverges : )3(r) ^ oo in a finite time. The non-perturbative solutions, on the contrary, includes effects 
of special relativity, notably the fact that j3 cannot be greater than 1. It has the major effect to delay the moment of 
"coalescence" of the two body. Finally we plot the motion of one of the companion over a large range of /3, and 
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FIGURE 3. (left) The ratio of the delay-time over the dynamical damping timescale as a function of jS at first order. These two 
timescales become comparable in the relativistic regime. (right) The behavior of the velocity in time, using the first order solution, 
the second order one and the PN solution one. 

compare the results of the PN approximation, the first and second order approximations. This is shown in Fig. |4] 
Although the behavior of the PN solution is quite bad near j3 = 1, we see that the three solutions agree very well until 
a separation of order 4m (a radius 2m), which is not too much compared to the Schwarzschild radius. This might be 
seen as a good news for PN expansions. 
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FIGURE 4. (left) The motion of one companion, from /3 = .18toj3 = .91, using the first order result /i. (right) The ultrarela- 
tivistic phase and the comparison of the first order solution (external line), the second order one (middle line) and the PN solution 
(interior line). 



CONCLUSION AND PERSPECTIVES 



The aim of this work was to extract as much as possible analytical results concerning the two-body problem in 
Nordstrom's theory, which is without doubts the purely scalar theory of gravity closest to GR. Once again, these 
analytical results are not a priori of direct interest as far as the construction of relevant templates is concerned, but are 
interesting since they enable to estimate the validity of methods used in the two-body problem in GR. For instance, 
analytical results derived here should be useful to test the efficiency of numerical codes. Those analytical results are 
already summarized in the introduction. 

Of course, an exact solution of the equation of motion may be one day derived, and that will be of great interest. 
Further work on the validity of the QE-scheme may also be done, see Q\. It could also be interesting to examine the 
behavior of the PN expansion in the plunge phase, compared to a numerical solution of the full equation of motion. 
We could therefore check if the PN predictions are very bad ones in the ultra relativistic limit (as expected, since the 
small parameter goes to 1), or if, due to some non-trivial cancellations (note for example the alternate signs in Eq. ( fT9] l 
and Eq. (|20]|). the PN picture reveals itself to be a good one. Such an alternation has already been observed in PN 
expansion of GR, see 1I21I1 . This behavior has not been explained yet, and we suggest that it might be understood in 
Nortsrom's theory. 
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